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We calculate the decoherence time of the ground state wave function of a nucleus in a high energy 
heavy ion collision. We define this time as the decay time of the ratio TrD^/(TrD)^ of traces of the 
density matrix D. We find that this time is smaller or equal to 1/Qs, where the saturation scale 
Qs is defined within the color glass condensate model of parton saturation. Our result supports 
the notion that the extremely rapid entropy production deduced for the early stage of heavy ion 
collisions at collider energies is to a large extent caused by the decoherence of the initial-state wave 
functions. 

PACS numbers: 25.75.-q,13.85.-t 


I. I. INTRODUCTION 

The physics program of the Relativistic Heavy Ion Collider (RHIC) has produced many intriguing results and posed 
a number of unexplained problems. One central question that has emerges is: How can hydrodynamical behaviour, 
implying local thermal equilibration and complete decoherence of the initial state, occur on a time scale which is 
considerably shorter than 1 fm/c? HQ. 

In an earlier article H we showed that the entropy per rapidity interval produced by decoherence alone is propor¬ 
tional to {RQsY, where R is the nuclear radius and Qs is the gluon saturation scale The resulting entropy per 
unit rapidity interval is of the order 1000—2000, which amount to a substantial fraction of the total produced entropy. 
We than argued that the decoherence time has to be of order l/Qs, as this is the natural scale of the process. The 
purpose of the present article is to substantiate this claim by means of a quantitative calculation. 

The fact that entropy produced by decoherence can play an important role in high-energy heavy-ion collisions was 
to the best of our knowledge first discussed in |^. 

Our article is organized as follows: In Section H we explain how we attack the problem and how we define the 
decoherence time. In Section HI we describe the calculation of the density matrix for gluons, which undergo a hard 
scattering process, in detail. Our calculation is based on work by Kovchegov and Mueller [iflIIS gluon production 
in heavy ion collisions. We generalize their derivation of probabilities to the level of density matrices. In Section IV 
we calculate the decoherence time using the results from Section HI. 


II. STRATEGY 

Our goal is to obtain an estimate of the decoherence time of the gluon distribution in a large nucleus (1), when it is 
hit by another very energetic large nucleus (2). In principle, the decoherence process is encoded in the time evolution 
of the density matrix I? in a very simple manner: For vanishing off-diagonal matrix elements the system is completely 
decoherent. In our case the dominant degrees of freedom are the gluons and thus the relevant density matrix is that 
of the gluons in nucleus 1. Decoherence thus manifests itself in a gradual disappearance of its off-diagonal elements 
Dki,k{ with the momenta ki ^ k'^. 

This sounds simple enough, but actually calculating the time evolution of off-diagonal elements in a complex multi¬ 
particle state is an extremely difficult task HQ. Therefore, one has studied so far mainly very simple toy models, 
like one harmonic oscillator in a termal bath of other harmonic oscillators, or rather specific situations, like neutrino 
oscillations Q. 

Luckily, also the situation encountered in high energy heavy ion collisions is such a special case, for the following 
two reasons: 

i) All gluons which undergo scattering are boosted into a part of phase space which was originally empty. 

ii) The description simplifies strongly in the rest system of one of the colliding nuclei. For an arbitrary Lorentz frame 
we do not know how to model the degree of coherence before the collision in both nuclei. In the rest frame of one of 
the nuclei, however, both nuclei can be approximated by different asymptotic descriptions. For a very fast moving 
nucleus one observes saturation and can describe the gluon field correlations along the lines of HI mi , while for the 
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FIG. 1: The chosen frame and notation 


nucleus at rest we can assume a simple isotropic Gaussian correlation of gluons with a typical virtuality of order 1/A. 
The precise form of this correlation does not really matter, because we will only study the correlations in that part 
of phase space mentionned in i), which is originally empty. 

Therefore, the elementary process of hard gluon scattering is sufficent to estimate the decoherence time. 

We perform our calculation in the rest frame of nucleus 1, which is hit by the highly Lorentz contracted nucleus 2, 
which is moving to the left. At lowest order in the strong coupling as the time evolution of the gluon density matrix 
is governed by the reaction Gi + G 2 ^ G'l , where Gi denotes a component of the gluon wave function of nucleus i 
(see Figure P). 

We investigate the gluon density matrix D{t) in a plane wave basis. Because the ground state wave function of 
nucleus 1 is a fully coherent bound state, the initial density matrix D{0) is strongly non-diagonal in this basis. As usual 
for a fast moving projectile, it is convenient to use light-cone coordinates {x+,x-,xj_) with x^ = (a;° pmx^)/y^2). 
Since, in our convention, nucleus 2 moves to the left, = {k^ — k^)/^/{2) = k 2 + is large and the nucleus remains 
localized in x~^ = x_, which thus denotes the light-cone position of nucleus 2. On the other hand, x~ = x+ « 
takes on the role of the “light-cone time”, by which the progress of the collision is monitored. This is also the starting 
idea of the Color Glass Condensate (CGC) approach |^, which is based on the assumption that most of the gluons 
in nucleus 2 hardly evolve on the time scale of the collision, endowing the gluon distribution with a “glassy” nature. 

The “time” evolution of D due to the three-gluon interaction is then given by 

D{dx+) = U{dx+)D{0)U'<{dx+) 


with 


U(dx+) = exp(—ii7intdx+/-\/2). 


( 2 ) 


For what follows, the variable r stands for the multiple quantum numbers completely specifying a gluonic state in 
nucleus 1: r = {fc,e^,6}, with wave vector fc, polarization vector e^, and color index h. The matrix elements of TJmt 
between two plane wave gluon states of nucleus 1 are then given by 


H, 




{ri\Hint\r[) = [ d^x±dx- [ gi(fc 2 +x--fc 2 ^-xr)ei[(fci+-fcj+)x--(/ci_L-/cG)-^v] 

J J 

9 f 0 .bc 9ct0 (^2 ^ 1 )-^ 901 (^1 ^l)a ^ 2 )^ 

[2fci+fcl+F2]-i/2gdg/72Tr [T‘^Alik2±,k2+)] ■ 


(3) 


Here 1/y v^^n-F is the normalization factor for a plane wave gluon in the light-cone formalism. V is some appropriate 
(—,T) normalization volume, which we will specify later. 
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describes the interaction with the incident gluons contained in nucleus 2, which moves very close to the 
speed of light. Since these gluons originate in many different nucleons, it is a good approximation to consider them as 
uncorrelated and to describe them as a Gaussian random ensemble representing an incident stream of Weizsacker- 
Williams (WW) gluons, all of which have positive /c 2 + momentum. 

The gluon fields of nucleus 1 will be characterized by a Gaussian ensemble of off-shell gluon fields with a spatial 
coherence length A of the order of the confinement scale Aq^^: 

= CabBxp {-[{x'_^ - x±f + (x'_ - x-f + (x'^ - x+f‘]/X^] . (4) 

We emphasize that we do not assume the gluon field in any nucleus to be dilute. We only assume that the typical 
momentum transfere fc_L is sufficiently large so that we can treat the interaction perturbatively in the coupling constant 
as=g^/47r. 

The WW gluon fields are given in eqs. (1-3). For our calculation it is crucially important that the color charge 
densities p°‘{xx,x-) and path ordered factors 

. J>e.p - xxl rt) (5) 

are uncorrelated for different values of x_: 

{p°‘{xx,x-)p'"{zx,z-)) = pAr(x±, x_)(5(x_ - - za_). (6) 

Zly qTT 

This and eq. (48) of [l^ allows to perform the calculation. 

Let us stress again that the asymptotic descriptions 0 and eq. (48) of treat the two nuclei in asymmetrically. 
Therefore, one cannot expect the final expressions to be manifestly boost invariant. Our calculation of the decoherence 
is specific to our selected Lorentz frame. While a manifestly boost invariant treatment of the decoherence process 
would be desirable, such a treatment would have to rely on a Lorentz covariant representation of the nuclear ground 
state. We do not address this interesting problem here. 

Because the incident nucleus 2 is a color singlet, its glue field vanishes on average: (x4“ (/c 2 )) = 0. This implies the 
absence of a contribution to the time evolution of the diagonal elements of the density matrix D in first order in g. 
The first nonvanishing term is thus of order corresponding to the term proportional to (x+)^ in (^, which is of 
second order in the WW fields of nucleus 2. 

The leading term for the evolution of the density matrix arises from the last term in eq. O: 

DrAr) = DrA0) + Y.'^-x'A)Dr>A0)WAA). (7) 

r' ,r 


where 



dx-\- (^-t-) ■ 


( 8 ) 


Our aim is to calculate the density matrix in the region of phase space populated by the outgoing gluons (momentum 
k{l)). This density matrix contains all the crucial information: 

• The ratio 


Tr {DA 
(Tr {D}y 


(9) 


is a measure for the coherence. For a pure state it is one, for a completely decoherent state it is much smaller 
than one, of the order of 1/N for a region of phase space with N states. We shall show that in the transverse 
directions (i.e. except for fc+) we get nearly complete decoherence. We will also see that this ratio depends on 
the observation time and we shall define as decoherence time Tdeco as the time after which it reaches 1/e. 

• The entropy of the final state can be calculated from 


SAf) =Ay {D{Tf)\ogD{Tf)}. (10) 

We argued in 0 that this is a sizeable fraction of the total entropy produced in a heavy ion collision. It should 
be possible to evaluate Eq. |Tn|l for the density matrix we obtain, but the calculation is highly non-trivial. 
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Therefore, we leave this task for a future publication. Having a model for the density matrix in hand it should 
actually also be possible to calculate other quantities specifying the initial state for the time evolution of the 
high-temperature phase produced in such collisions. 


The symbol (• • •) indicates and average over the WW fields of nucleus 2. 


III. THE DENSITY MATRIX RESULTING FROM HARD GLUONIC INTERACTIONS 

For the WW fields of nucleus 2 we have ^ 2 - = 0, because the fields A{x) do not depend on x+. We therefore get 


(A“ (fc')A“(fc2)) = / 


r4 / ~^ — ik'oX' -\-lk 2 X j \Ol ( f\ 




( 11 ) 


= (27r)2<5(fc^_)5(ft2-) J 
The relevant term is now: 

Wii.i'i = 

= J dk'^_ J dk2-5{k'^_ + k[_ - ki-)S{k2- + ki- - k-) 

J dk'2+d^k'2xeik'2+) J dk2+d^k2±0ik2+)Sik'2+ + k[+ - ki+)S^{k'2^ + k[^ - kix) 

S{k2+ + ki+ - ki+)S‘^{k2± + ki± - ki±) 

gfabc fffa/3 {k 2 - k[)^ + {k[ + kl)^ + g-ya {-ki - 


all 


abc 


9&P (fe - ^i) _ + a^j (ki + ^i) - + 9i& (-fci - fe) 


€ 

T’^Al{k' 2 ^,k' 2 +)] 2Tr [rMf fe+)l. 


2Tr 


( 12 ) 


Here the first line (the k 2 - integrals serves more or less only as reminder.) Note that there are no normalization 
factors for the incomming gluon states, i.e. no factor [2fc(_|_fci+U^]“^/^. The reason is that our initial states are 
highly virtual, interfering gluon states, i.e. no momentum eigenstates. For the time being we absorb all normalization 
factors into our definition of d{x,x') which therefore is not dimensionles but in coordinate space has the dimensions 
Energy^. We do not know how to properly normalize a density matrix for an arbitrary virtual field configuration. 
Luckily we will only need the normalization of the diagonal elements. 

We first focus on the ensemble average of the expression in the last line, involving the WW fields in nucleus 2: 


(J) = U 


2Tr 


A:' + )J 2Tr (fc 2 ±, fc 2 +)J 

{wj_ - y_L)“ (w_L - yx)°‘ 


d^y±dy-0{y- -W-) d‘^y±dy-9{y- - W-)- 


|w_L - \w± - y_Lp 


2Tr 


T'^Soiwx,y-)T<^So\wx,y-) 2Tr T-So{wx,y-)T'^So\wx,y-) /(y±, y-)/(yi, y.)) . (13) 
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We now assume complete factorization of the average, i.e. 


'2Tr 
= /2Tr 


= (2Tr 
a, 


X ■ 


'2ttN, 

Using the identity 


T‘^So{w±,y-)T‘^So\w±,y-)\2TT[T^So{w±,y-)T‘^So\w±,y-)\p‘^{y±,y-)p‘^{y±,y-) 
T-So{w^,y-)T^So\w^,y.)]2TT[T-So{w^,y-)T^So\w±,y-)\){p^{y±,y-)p^{y±,y-)) 
T^^Soiws,, y.)T‘^So\w^ , y_)] 2Tr [t^Sq i^v^ , y.)T‘^SoHw±, y-) 

PN{y±, y-)S^‘^S{y- - y-)6^{y± -yi_). 


2Tr 


rjiCi j^erp€ 


2Tr 


= = 2Tr 




to rewrite the product of two color traces as a single trace, we get 

{J) = U-2 [ / d^w±dw-e'^'^^+^--'^^^-^^^S{k'2_)6Ck2-) 


(14) 


(15) 


d'^yj_dy-9{y- - W-)9{y_ - w-) 
2Tr 


{w± - yi)°‘ (w L - yA Y 


\w±-y±Y \w±-V±\'^ 2 'kN, 
SY\w^,y.)T-S^{w^,y.)SY\w^, y-)T^So{w^,y-) 


■PN{y±,y-) 


(16) 


Now we use eq. (47) from na and assume that for a ^ d the ensemble average vanishes. The fact that the original 
expression (47) has Sq and factors interchanged and thus is the complex conjugate of our expression, does not 
matter, because the result is real. 


Tr 

gaa 


SY\w^,y.)T-Soiw±,y.)SY\w^,y.)T-Soiw^,y.) 

]s^2_i SY^{w±,y-)T°-So{w±,y-)SY^{w±,y-)T°-So{w±,y-) 'j 


(5““C7fW 

m - 1 


■ exp -g 


I "^Prel-^c 
- 1) 


a;G'(a;, |rZ;_L - u;_l| ^)(w± - w±)^(y-+ y®) 


(17) 


where B? — 

This gives the following result: 

(J) = U-2 f / d^w^dw-Y^'^^+^--'^^^-^^'^S{k'^_)6Ck2-) 


Jd^y. 

2CfN, 

m - 1 


( 0 ) 

rv- 


'-y 


( 0 ) 


dy-9{y- - w-)9{y_ - w-) 


{w± - y±)°‘ {w± - y±)°‘ as 
|w_L - y±Y l^ej- - yj-P 27riV, 


exp -g 


2 Prel^c 

^ 4(iV2 - 1) 


xG{x,\w±-wj_\ ^)(w_L - W±)^(y-+y®) 


■PN{y±,y-)S° 


Next we perform the y- integration. To do so we make one more approximation. I assume <C W-,W- and 
substitute the ^-functions by 9{—W-)9{—w-). 


(18) 


iJ) = 


C /*0 r pQ 

y-^ / d^wi / / .72 


d y±S 


d w± 

^ —OO 

a(w_L -y±)“ (w± -y±y 


d'^w± I dw-ed^^+^—’^^^-^^'>5{k'^_)5{k2-) 

J —OO 

4(N^ - 1)27/® 


\w±-y±\'^ \w± - yFl"^ g‘^TrNcPrei\w± - w±\‘^xG(x,\w± - w±\ 2) 


1 - exp -g- 


I 2J/P VpreWc 
4(iV2 - 1) 


{w± — W±YxG{x,\uI± — w±\ ^) 


2TrN, 


■PN{y±,y-) 


(19) 
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In the spirit of the discussion above we assume 

PN{y±,y-) 

2j/i°Vrel 

and use the relation (valid at leading logarithmic accuracy) 
dxG{x, as{Nc - 1) 


= 1 


91nQ^ 


2ttNc 


xG{x, Q ) = 


,2\ _ 1 ) 1/^2 / ,,2\ 


2TTNr 


■ln(Q //i ) 


to obtain 


(J) = -1/-2 / / d^y_e-Hfc2+™--fc2^-»-L) / d^w± / 


d y± 


-4 


iw±-y±)°‘ _ 

\w±-y±\'^ \w±-y±\'^ g^ 7 TNciw± - w±)'^ln{p,‘^\w± - w±\‘^) 


2 


Next, using the two dimensional Green’s function 


{w± — w±yxG{x^\w± — w_}_\ ) 


5{k'^_)5{h_) 


{w^ -ys.T . f dh± 

|w±-?/±|2 J 2tt q\ 


the y_L integral can be performed: 


d^yx- 


(w± - y±)^ (w± - y±)^ _ _ f ^ 
\wx-y±\^ \wx-y±\^ - 


2tt 2tt 




— f d^ai gig-L-l^’-L-^’-L) d±Q± 
-Jdq^e 


dq± 


= / n± I d(t)e 

Jo Q-L Jo 

= 2.S--rdgx 


igj_ |u;_L —-mx I cos 0 j COS (j) 

sin^ ch 


_ J{(g±|w± - w±|) 

10 q± ’ V'^i('I-Lk-L - +-^2(9^1^^ - W_l|) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


Up to the logarithmic divergence at gj. —> 0 this integral is independent of \w± — as one can see by substituting 
q±\w± — —> z. To regularize the IR divergence we introduce a lower integration boundary for z in the form 

p,\w± — and insert the finite z ^ 0 limit of the Bessel functions: Jo(0) = 1, J 2 {—) = 0, obtaining: 


/ fl\W± —Wd_ 


Z 2 i 


(25) 


We thus finally end up with the result 


(J) = U-2 / d^w^ / / dJw^ / rf^Zi_eHfe2+»--fc2^ *_L)^aa^oa 


g'^Nc\w^ - w_l| 


^ ( 2‘^y-^'^P^e\Nc \2 r<( |~ I-2n^ 

1 - exp y-g 4 (jy 2 _ 1 ) W-L “ "^-l) 2 ;G(x, |u;_l - u;_l| ) 

I insert a suitable e-prescription to perform the integration over W-: 

25““(5““ f d^wxd^w± j(fc, , ,-fc,, .^ 5 , ) 


5{k'2_)5{k2-) 


(26) 


{J) = 


g^N,V^ J |wi-u.i|2 


pO pO 

' dw- / rf^;)_e-i((fc 2 +-eie)«;--(fc 2 +-ie)*-) 

— 00 J —00 


1 - exp -g 


2 2yi° VpreWc 

4(NJ - 1) 


('W± — W±)^xG(x, jlCj, — W_L 


-2n 


dikUdik-) 


25““5““ 

g^N,V^ 


d^w^ — 


/ d'^w^- " 


■^2+ 


le 


■ k2+ - le 


-e |w_L — ^ 


1 f 2 ^ 2 /— '^PrelNc , ^ ^2 | ~ |— 2 \ 

1 - exp I -g {Wj_ - w_l) xG{x, |w_l - w±\ ) 


5{kk)5{k2-) 


(27) 
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Next we substitute Zj_ = w± + w± and z± = w± — w±, giving 


iJ) = 


(A;'_^ + ie)(fc2+ - ie) J 
2 2yi° Vprel-ZVcZ^ 


1 - exp -g^ 


m! - 1) 


' (fz±j d2^^gi[fe2i'(^±+z±)/2-fc2J_'(Zj_-Zi)/2] 

xG{x,\z±\~‘^) ] S{k 2 _)d{k 2 -) 


(28) 


Now we can insert the definition of the saturation scale from eq. (17) of m In doing so we identify r in that equation 


with : 


This gives: 


. , ,_2^ Qlzl 

4 (iV 2 _i) ) = ^ 


{J) = 


2^-^-(27r)^ 6^{k!2^-k2±) 

g^NcV^ (^2+ + ie)(fc2+ — ie) J 


1 — exp — 


Qlzl 


6{k'2_)6Ck2-) 


(29) 


(30) 


To do the z_L-integral we use the fact that the integral is well behaved at —*■ 0 and we assume some e prescription 

to make it convergent at ^ oo. 


_ f <Pz±_ -u 


K= I i_exp(- 


Substituting z_l ^ z_l + i2fc2_L/(QsU), we get: 

-,2 ^1 r 


Qlzl 


Q 


2 ri 


= ^ / du d^z±e 


<2 ik2j_-z± 


exp - 


Qlzlu 


K = 


Qi 


du / d z_L exp — 


Qlz^u _ k'^x 

4 Q2y 


du 


k' 


• exp 


2± 


Qlu 


(31) 


(32) 


Next we substitute u k' 2 ±/{Qlt) to obtain: 

K = n 

yielding finally 


dt t / k'o I 

I —e-* = ttEi ' 2J- 


Qi 


{J) = 


d°°d“°(27r)3 5'^{k'2^-k2±) 

g'^NcV'^ (^2+ + ie)(fc2+ - ie) 


El 


^ 2 ± 

Qi 


SikUSih-) 


(33) 


(34) 
























thus we get: 


Wij i,j = J dk2_ J dk2-6{k2_ + k[_ — ki-)6{k2-+ ki-— ki_)S{k2_)S(k2-) 

(2^)^ ^ 1 ^. 


V 


-fabcfabc / ^^ 2 + d?k’22_e{k'2+) 


dk2+d^k2A_9{k2+)5{k2^ + k[^ - ki+) 


2ki+ki+ 

<^^(^ 2.1 + K± - ki±)d{k2+ + h+ - ki+)S^{k2± + h± - ki±)6‘^{k2^ - h±) 

9af} {k2 - k[)^ + gfs^ {k[ + fci)^ + (-fci - 

5/ (^2 - fci) _ + 9p;y (ki + fci) + 57 “ (-^1 - ^ 2 ) 


q;=1,2 


(fc^+ +ie)(fc2+ -ie) 


El ( ^ ) e(fc;A*(fci)^(^e(fcire*(fcir 


{2'K)^S^{k[^ - ki± + ki± - ki±)fabcfabc 


V 


S(k[_ — fci_)(5(fci_ — ki-) 


I dfc' + 0(fc' + ) / dk2+0ik2+)^== 

\j2ki+ki+ 

9af3 (^2 - ^ 1)7 + 5/37 iK + ^l)a + 57« (“^1 “ ^ 2 )/? 
5 / (^2 - fci) _ + 9pj {ki + fci) + 57 “ (-fci - fe) 


= (5(^2+ + ^1+ - fcl+)5(fc2+ + fci+ - fci+) 


a=l,2 


-El 


l/2 ' 

^ 2 L 


(^ 2 +-l-ie)(fc 2 +— ie) Qs 
( 27 r) 3 ( 52 (fc'^ - ki± + ki± - h±)fabcfabc 


{k[)^e*Ckir(j2<k,y^*{kr 

e 

0 (fci+ - fci+) 0 (fci+ - k[_^)S{k[_ - fci_)5(fci_ - fci_) 


Vy2ki+ki+ 

3 a /3 {^2 ^1)7 “I" 5/37 (^1 “I” ^i)q; 57a ( ki ^2)/? 

5 / (^2 - fcl) . + 5/37 (^1 + ^0 + ^ 7 “ (“^1 “ ^ 2 ) 


El ( ^ ) e(fci)^e*(fci)^(^e(A:ire*(fcir) 


(fci+ - k[_^_ + ie){ki+ - ki+ - ie) 


Q 


(35) 


with k 2 j_ = ki± — k[j_ and k 2 ± = ki± — ki±. 

The WW-gluons are on-shell, i.e. k 2 _ = fe- = 0 and thus fci_ = fc(_ and fci_ = fci_. Therefore we can substitute 
in general /c^ = fci — k'l and k 2 = ki — ki. 


W 


ii,i'i 


(27r)3j2(fc/^ - k, + ~k,^)f,i,cf,is 

Vy2ki+ki+ 

gap {ki - 2k[)^ + 3/37 (fc( -b fci)„ -b 37 a (- 2*1 -b *()^ 


9 ^a {kl ~ 2 fcl^ ^ + 5/37 (~ 2 fci -b *1^ . 


(* 1 + - *(+ -b ie)(fci+ - * 1 + - ie) 


-El 


k'l^ 

Ql 


,{k[fe*{hf{^e[k,rc*{kr 


(36) 


At this point we have to make an assumption about how far in rapidity a gluon of nucleus 1 is scattered. Let us 
require an average rapidity T > 1 of the scattered gluons (see Eig.|2l. Decoherence will be effective, when Y is large 
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enough that the final phase space is originally empty or only sparsely occupied, which is the case for Y > 1. This 
condition also implies /ci+, fci+ ^ k[^, fci+. Decoherence may also occur inside the target region (Y < 1), but we are 
not concerned with this question here. 

With these assumptions it is now possible to greatly simplify Eq. In the square brackets we keep only ki^, 

ki o and fci.s, ki^^. We neglect ki^a, because a denotes a transverse direction. We can also drop ki^^ and ki^^y because 
of the projector at the end of Eq. Thus we are left with 

ki ■ ki 

= 8 fci+fci+e(fcj)_ e*{ki)- ( 1 + l 

V 2fci • ki I 


= 8 ki ■ e{k[) ki ■ e*(fci) 1 + 


feyj. ■ fci,± 

2ki ■ ki 


Making use of: 


^e(fcire*(fcir 



kfkj \ 

ki -ki) 


The next task is to evaluate the correlator 


(38) 


(e(fcl)_ (39) 

Next we have to evaluate the convolution with the density matrix of the product of polarization vectors and the color 
factor: 


e (^1 ) —)polarization,color ^ EE e(fcl)_ e*Cki)-fabcf,i,D(^b,b,e{k[),eCk[);k[,ki) (40) 

6,6 

To do so we Fourier transform into coordinate space and simply assume that there (in the rest frame of nucleus 1) 
the correlation of the polarization vectors and colors is given by a simple Gaussian with correlation length A, i.e. we 
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assume that I? is a Gaussian in coordinate space. 

ifabcfabc^iK)- e*(fci)- ) = CScc J (fy± d^y± dy'_ dy- dy+ dy+ 

exp |i(/ci_^ • yj_ - h,± ■ y± - K^+V- + ^+V-) “ -2/+ + ^-y+)} 

exp I - [{y'^ - yif + {y'_ - y_f + (y^ - y+f]/X^^ 

= J d‘^Y^+^ d^Y^-^ - exp { - ^[{k[ - ki) ■ F(+) + {k[ + ki) ■ 
exp{[(F|-V + + (f[-V]/A"} 

= C6,&{V^X)^6^{kl^ - ki,±)S{k[_ - k-)S{k[^ - ki+){27Tk 

exp I — + {k'i_ + fci-)^ + {k[_^^ + fci+)^]/4| 

= C6,,iV^XkSkK± - ~ki,±)S{k[_ - k-)diK+ - k+){2Trk 

exp I - X'^[{k[^j_f + {k[_f + (fci+)^]| (41) 

with F(+) = y' + y and y(“) = y' — y, C is a constant to be determine from a suitable normalization condition. 


Now we are at the point where the normalization of d, respectively D, has to be discussed. To do so we start from 
the energy associated with the quadratic part of the Yang-Mills Lagrangian, 

C = -^id^A:-d,A;)d^Aa., (42) 

which should be equivalent to some integral of the form 

B = /(^ -'i(p.p) («) 

with a one-particle, not yet normalized, density matrix d{p,p). ^From Eq. 14211 we get 

E=^Jd^x (^id°A^kx)f - id^A^kx)f + id^A^kx)){d^A^kk) - (5M*'^(x))(a^A4^*(x))) (44) 


Substituting 


A^f^{x)A^''{x) 


d'^p f d'^p 


( 27 r )4 J {2 ttY- 


dp%b,b-p-p) 


(45) 


we get 


E = 


1 f d^p dp^ 


2 J (27r)5 


(^(pV + ip)‘^)d^kb,b;p°,p;p°,p) - {p)‘^d°°{b,b]p°,p]p°,p) - p’'fd^Yb,b;p°,p;p°,p)^ (46) 


For simplicity we choose the gauge p^A^ = 0, allowing the last term to be re-expressed in the form p^k^dP^. Our 
model assumptions Eq. ED imply that we only need the diagonal matrix elements of H. To fix it we require that 
Eq. 14611 reproduces p° for a plane wave gluon with 4 momentum p^, i.e. we impose: 


f = 


d^p dfP 
(27r)5 


,p°p° 


{p) )dp{b, b;p°,p,p°,p) - (pV + (p) )df{b, b;p°,p,0 


\P)) 


(47) 


for a density matrix of the form 

dp‘'{b,b-,p°,p,p°,p) = 5^k5vk{2TTf5'^(j)-p)f{p°,p°,p) 


The solution for / is obviously 


P^Att 


(48) 


/(p0,p°,pl=^(p0-p°) 


(p0)2 + 


(49) 
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This motivates our assumption that the diagonal matrix elements of our d are related to those of the properly 
normalized density matrix D by the factor (setting (5(0) = T) 


df^’^{b,b;p^,p,p^,p)=T 




Df^’'ib,b;p^,p,p^,p) VT, 


(p0)2 + 

where D is according to Eq. 6U in our model given by 

D^uib,b]p°,p,p°,p} = -^TT^A'^expl - X'^[{p±f + {p-f + {p+f]^]^5^,u\^J.,v^±■ 


(50) 


(51) 


Here the superscript “_L” indicates that the Kronecker symbol contributes only for the transverse directions p,iy = 1, 2, 
and the factor A encodes the unpolarized nature of the gluons in the target nucleus. The normalization of the trace 
of the density matrix demands that 


VT 


d'^p 


f,{b,b-,p°,p,p°,p) = 1 


(52) 


for a one-gluon state. 

For a many-gluon state this should be normalized to the total number of gluons which is an ill-defined quantity. 
We therefore choose to substitute for this general case E = hy the total energy of the gluons Ec{Ql) in nucleus 1, 
at the transverse scale with which this nucleus is resolved, which is the saturation scale. 

Combining everything the density matrix in the boosted, previously unpopulated region of phase space is 


(C) c', ki, ki, , ki) — 


d% f d% 


(27r)4 J (27r)4 


{2TryS'^{k[ — ki) 


VT 


Wii.ii 


47rE;G(Qf)T 

(fc?)2 + (fci)2 


D^,{b,b;k[,h) VT 


(53) 


with 

fabcfabceiK)- fabcfabc ^ NJcc^ . (54) 

The factor | is motivated by the assumption that in a nucleus at rest (the target nucleus 1) all gluons are so 
highly virtual that the transverse and longitudinal polarization components contribute equally. The gauge degrees of 
freedom, of course, do not contribute. 

Combining all our results we get for the one particle density matrix in the final state phase space region (substituting 
ki by p for notational simplicity): 


D, , = 


f d'^p iTTEGiQs)l6n^X'^NJcc ( ^2^^ 2 , 2 , 2 11 

J w w^{{p«y‘ + m «p I - ^ + p+1} 

(27r)3(5^(A:i_L - ki±) 


0{ki+ - p+) 0{ki+ - p+)S{p- - ki_)5(p- - fci_) 


2fci+/c 


1 + 


8 I 1 -|- 


ki,± ■ fci,-L \ _1_ 

2fci • fci ) (/ci+- p+-I-ie)(A:i+- p+- ie) 


-El 


(^i± -Pi.f 

Ql 


(55) 


In order to evaluate this expression, we note that we are interested in the scattering of gluons into states with 
fci+, ki+ ^ ~ 0{X~^). We can then neglect the p+ dependence of the two denominators, drop the step functions, 
and obtain: 


IGtt^ 


e-^ iP-L+P-+pV 


(27r)4 pI 


S{p_ - fci_)d(p_ - fci_) El 


(fci L -p±)" 


Ql 


= ^(5(fci_ - 


poo nOO pOO / 

fci_) / di^ / dp+ / ( 

J J— oc J—oc \ 


{ki± -pj)^ 


5{ki.-k-) J 

Q2 / 


7r3/2 

2 ^ r°° 

= -j=5{ki- — fci_) I d^ e 


/TT 


L I 


exp 


C2u2 




(56) 
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We now have traded 3 integrals (over and p±) for two integrals (over ^ and u). This may not seem like much 
progress, but it turns out that the integral over u can be done after the substitution s = + 1): 


L I 


exp 




21,2 
1± 


+1 


ds 


3 ^^1±/Qe 


I s2 ^2Q2 
qW^ 


1 


El 


ki^e 


Qie V \Qie + ^ 


-Eiiki^e) 


We finally substitute ^ ^ in the remaining integration and obtain for the expression 


(57) 


^1,1 = 


(47r)4A3£;GiVAc 5\ki^ - ki ^)5{ki. - h-) A ^ ki^j_ ■ 


3^2 


2fci+fc 


1 + 


2 /’°° -\^k^ 

' T^k 


El 






2ki ■ ki 

-Eiix^kf^e) 


(58) 


This expression describes the density matrix of the liberated gluons, which are scattered out of the target nucleus 
1 by the quasi-real gluons of the fast moving projectile nucleus 2. We note that the density matrix is diagonal in 
the momentum components fci_ and ki±, but not in the component fci+. The physical reason for this asymmetric 
behavior is that the projectile nucleus is moving very fast in the x_ direction. This implies that the distribution of 
its gluons in is very broad and leads to interference of excitation amplitudes of gluons from the target nucleus into 
final states with different values of ki+. 


IV. THE DECOHERENCE TIME 


We now calculate the ratio for the density matrix from Eq. (CTl . we define 


F(fci_, fcij.) = 


e 


x-^ki^e 




(59) 


Because any constant factors will drop out of the ratio ®, it is sufficient to consider the fc-dependent part of j, 
which we call the reduced density matrix: 


I?i,i = _ ki.)5^{ki^ - ki±), 


(60) 

'ki+ki+ 

where we have neglected the factor (1 -|- ...) deriving from the polarization sum, which is of order unity. We obtain: 


TrX> = 


VT 

VT 


/ 77 poo poo 

dki.J d^ki^F{ki.,ki^) 


S{0.)6^{0j_) 


dkl+ 2 


ki+ 4A 




(61) 


The details of the integration can be found in appendix A. For TrD^ we obtain 

2 „ „ 


TrD^ = 


VT 


d^k / d'^k- 


1 


ki+ki+ 

S2 


F{ki_,k±)S{ki- - ki_)S^{ki± - ki±)F{ki_, k±)S{ki- - ki_)S^{ki± - ki±) 

& I '-t I 


Which can be simplyfied using the identity (see Appendix A) 


f dkfj_Ei{aklj^)Ei{bkfjJ = -\n —-1—In 

J —oo Odd 


( 63 ) 
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FIG. 3: The function I(z) and the ratio I{z)lz 


which leads to 


TrX>^ = 


( VT \ r dki+ r dki+ 2 .^ f 

J T7i IT el 


d^ rdx 1 


i,j=l 


>'1 + 

1 1 I 

— In- - H-In —;—- 




ibj 


a* 


Oi 


with 


Oi = 


= 


A2Q2^2 + 1 


We now define the integral 

I{XQs) = 


X^Qlx^ + 1 

di rdx I 


02 = x^Qie 

b2 = X^Qlx^ 


eJr 


E (- 1 ) 


i+j 


1 ai + bj 1 Oi + bj 

— In- - H-In —^—- 


(64) 


(65) 


( 66 ) 


which approaches zero for XQs 0 and approaches for asymptotically large argument the limit /(oo) = |(^/2 — 
l)ln 2 « 0.3828. One can see the latter as follows: XQs is much larger than unity and also > 1. Therefore 
02 2 > oi ~ 1 and 62 ~ 1 and the contribution with i = j = \ dominates. In this limit, the square bracket 

assumes a value of about 2 In 2. Substituting this value into the integral ( 68 ), we obtain the asymptotic limit mentioned 
above. The exact form of I {XQs) is given in Fig. El 
The degree of decoherence is given by the ratio 

TrI?^(totaZ) Tr [I? (target) + , , 

(TrT>(tota/))^ (Tr[I?(target) +1?])^ 

where 'D{target) is the density matrix of that part of the phase space which contained the original target. We are 
interested in the situation that most gluons of nucleus 1 have already undergone a hard scattering, i.e. T must be 
large enough, typically of order > 0.1 fm/c. Then we can disregard ^{target). 

Thus we finally obtain for the ratio which characterizes the degree of decoherence after the hard gluon scattering 

Trl?^ 16X^I{XQs)Q^ 8I{XQs)X tt 1 

(^^^ " 2A 7r3QfJ2(o^)^(o_) = 


( 68 ) 
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With the “observation time” 5(0_) = T and ^^(Ox) = A^/ TT, which is the correct normalization for eq. because 

VT 

We now define the decoherence time as the value of T for which the ratio becomes equal to 1/e. This gives 


_ 8eI{XQ,) 1 
0FQ.A ' Q, 


(70) 


For realistic values (A = 0.3 fm, Qs = 1 GeV) the first factor is numerically of order unity (more precisely it is 1.865) 
for these values) and we can conclude 


rdeco ^ ^ ("^1) 

This is our main result. The decoherence time is of the order of 0.3-0.4 fm/c, and thus large enough to neglect 
T>{target) in Eq. ll^ . and drops with increasing saturation scale Qs, i.e. with increasing collision energy. However, 
Qs increases so slowly with s that the latter effect is rather marginal. We use A(5s=l-5 and think that the uncertainty 
in this value is not large. Nevertheless, we want to point out that even if one varied this product by a factor of 
two, the ratio I{XQs)/XQs would change only marginaly and (I71II stayed valid. (For XQs much smaller than one our 
description stops making sense.) 


V. CONCLUSIONS 

We have calculated the characteristic decoherence time in high energy heavy ion collisions due to gluon scattering. 
We find that this time is substantially shorter than 1 fm/c. This result furnishes the remaining logical link in our 
argument that decoherence alone can explain a substantial part of the entropy production during the earliest phase of 
a heavy ion collison. We note that the decohered partonic state of the colliding nuclei is not yet thermally equilibrated. 
Additional interactions among the decohered quanta, such as those invoked in the bottom-up scenario of equilibration 
|l2j| or modified versions of it mill are required to achieve full equilibration. 
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VI. APPENDIX A 


We here evaluate the integral appearing in eq. (60): 


I' = 


dki- 


r§e-> 

, Ji r 


El 




-Eiix^ki^e) 



poo 

/ d^k 

2A 

J — oc 


rd^ 

2A 

Ji ¥ 


b/s 2’ 

2A 


1± 


d( 




1 ei + ^ 


x^e 




where we used the relation 



dzEi{z) 



(72) 


(73) 
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We also derive the identity (62): 

cP'ki^Ei {akn_)Ei {hki±) 


TT 

2 


dk 


dt 


2^ / _g-ofcfit 


ds 


— bk': , s 


'0 


TT 

2 . 


dt 


1 t 

poo 


Jl ^ 

ds 1 
I s at + bs 

at V ot + bs 


TT f°° dt , , . , M 

= ^ y t2 +^) “ 

'■°^dy 

i/b y 

In 


TT 


■ ln(2/ + 1) 


TT 

TT 

2 


j - iln(y + 1) 

y + iy y 


i/b 


0 V a 


a \ 0 


(74) 
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